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We present a large grid of computed far- and mid-ultraviolet spectra (850 A to 2000 
A) of the integrated light from steady-state accretion disks in luminous cataclysmic vari- 
ables. The spectra are tabulated at 0.25 A intervals with an adopted FWHM resolution 
of 1.0 A, so they are suitable for use with observed spectra from a variety of modern 
space-borne observatories. Twenty-six different combinations of white dwarf mass M w d 
and mass accretion rate m are considered, and spectra are presented for six different 
disk inclinations i. The disk models are computed self-consistently in the plane-parallel 
approximation, assuming LTE and vertical hydrostatic equilibrium, by solving simulta- 
neously the radiative transfer, hydrostatic equilibrium, and energy balance equations. 
Irradiation from external sources is neglected. Local spectra of disk annuli are computed 
taking into account line transitions from elements 1-28 (H through Ni). Limb darkening 
as well as Doppler broadening and blending of lines are taken into account in computing 
the integrated disk spectra. The radiative properties of the models are discussed, includ- 
ing the dependence of ultraviolet fluxes and colors on M w( j, m, and i. The appearance 
of the disk spectra is illustrated, with regard to changes in the same three parameters. 
Finally, possible future improvements to the present models and spectra are discussed. 
The synthetic spectra are available as machine-readable ASCII files via ftp. 
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1. INTRODUCTION 

In most varieties of cataclysmic variable (CV) 
stars, the mass transfer between the Roche-lobe filling 
secondary star and the white dwarf primary star oc- 
curs via an accretion disk. It is clear from decades of 
observations and theoretical studies that this disk can 
dominate the ultraviolet and visible spectra of CVs, 
at least for the more luminous classes such as novalike 
variables and dwarf novae in outburst. Observations 
of CVs with the IUE observatory established the im- 
portance of the disk light in the mid-ultraviolet (mid- 
UV; roughly, ultraviolet wavelengths longward of H 
Lyman-a), and Hubble Space Telescope (HST) obser- 
vations using HSP, FOS, and GHRS have resulted in 
improved data for many key systems, with improved 
spectral and temporal resolution and enhanced signal- 
to-noise ratio. Observations with the ultraviolet spec- 
trometers (UVS) aboard the VOYAGER spacecraft, 
with the Hopkins Ultraviolet Telescope (HUT), and 
with the Space Shuttle-borne ORFEUS experiments, 
have shown the same dominance of the disk in lumi- 
nous CVs for the far ultraviolet {far-UV; from the 
Lyman edge to Lyman-a). The disk may remain 
dominant in less luminous CVs (mainly the dwarf no- 
vae in quiescence), or it may contribute significantly 
less than the white dwarf, especially if the disk is 
highly inclined to the observer's line of sight. In one 
class of CVs, the AM Her stars, the accretion flow is 
channeled by strong magnetic fields and a disk is not 
formed. 

To make the most of the large and increasing 
archive of high quality mid- and far-UV observations 
of disk-dominated, high-luminosity CVs, and to dis- 
tinguish reliably between white dwarf and disk spec- 
tra in less luminous cases, it is clearly necessary to 
have up-to-date models for the spectral energy dis- 
tribution and detailed line spectra of disks, for a va- 
riety of accretion rates, central white dwarf masses, 
and viewing angles (orbital inclinations). As obser- 
vations have been collected over the years, different 
techniques, involving different levels of approxima- 
tion, have been applied to this task, and a few small 
grids of models have been produced. For example, 
Williams and Ferguson 1982 described a small grid of 
LTE disk models, spanning the transition from opti- 
cally thin to optically thick in the continuum, in which 
the disk was assumed to have constant temperature 
in the vertical direction at each radius. Wade 1984 
discussed the UV and optical continua for CV accre- 



tion disks, comparing two grids of spectra, based on 
blackbody energy distributions and computed stellar 
fluxes from the Kurucz (1979) grid, la Dous (1989) 
introduced some absorption lines from metals into 
her computed ultraviolet spectra of disks, which were 
LTE models based on the Eddington approximation 
without any attempt to ensure energy balance, la 
Dous also summarizes prior work by other authors 
(see also la Dous 1994). 

Recent data obtained using the instruments listed 
above, with their higher S/N ratios and spectral reso- 
lution, demand more detailed models, computed with 
some attention given to the properties of disk spec- 
tra that may distinguish them from sums of stellar 
atmosphere fluxes. The spectra in Wade (1984), for 
example, are coarsely binned and entirely inadequate 
for the resolution achieved in HUT or HST spectra. 
Many studies of individual CVs in recent years have 
made use of improved detail in model spectra (e.g., 
Long et al. 1994; Knigge et al. 1997; Knigge et al. 
1998), but there has been a lack until now of a grid 
of mid- and far-UV spectra of disks, presented as a 
body. Such a grid can enable comparisons among ob- 
jects, allow the range of temperatures in the disk to be 
estimated, help identify lines and blends that can dis- 
tinguish disk spectra from white dwarf spectra, and 
perhaps in favorable cases help determine the mass 
accretion rate and orbital inclination in some CVs. 

It is also desirable to have a new grid that explic- 
itly tabulates spectra for disks viewed at a variety 
of inclinations, since it is clear that limb darkening 
of the UV light from a disk is of comparable impor- 
tance to purely geometric effects for optically thick 
disks viewed at high inclination (e.g. Diaz, Wade, 
& Hubcny 1996; Wade 1996; Robinson et al. 1995; 
Robinson, Wood, & Wade 1998). This was shown 
in the early calculations of Herter et al. (1979), and 
has been illustrated from time to time (e.g., la Dous 
1989), but most discussions of disk spectra since then 
have treated angle-averaged quantities. Diaz et al. 
(1996) discussed several aspects of limb darkening for 
a grid of 14 model disks, but tabulated limb-darkened 
fluxes at only one wavelength. 

In this contribution, we present model mid- and 
far-UV spectra for accretion disks in CVs, for a range 
of white dwarf mass M w( j, mass transfer rate to, and 
inclination angle i, that covers the expected range of 
interest for luminous (or "high-state") accretion disk 
systems. These models treat steady-state accretion 
disks in vertical hydrostatic equilibrium, with the en- 
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ergy balance of the disk computed in a self-consistent 
way. The modeling assumptions and techniques are 
described in Sections 2 and 3. Properties of the model 
disks and especially their spectra are the subject of 
Section 4. Such grids of model disk spectra need not 
be perfect to be useful, especially when used in a 
comparative way, or as a starting point in examining 
which way models should be improved or extended. 
We discuss some of these issues briefly in Section 5. 

The spectra discussed here are available in tabular 
form as ASCII files via anonymous ftp (see Section 
3). 

2. MODEL ASSUMPTIONS AND CHOICE 
OF GRID 

2.1. Steady-State Disk Model Computations 

Our grid of spectra is constructed as follows. The 
disk is assumed to be axisymmetric and geometrically 
thin. The disk is assumed to be in steady state, that 
is the mass transfer rate m is assumed to be the same 
at all radii. Each elementary area of the visible disk 
surface is assumed to make a contribution to the in- 
tegrated spectrum of the disk, according to its pro- 
jected surface area and at a projected Doppler shift, 
computed assuming circular Kepler orbits around a 
central white dwarf star of mass M w( j and radius i? w d- 

The effective temperature, which describes the to- 
tal radiation energy flux at the disk surface, is given 
by (see, e.g., Pringle 1981) 

T eS (r)=T*x- 3 '*(l-x- 1 / 2 ) 1 / 4 

where x = r/i? wd , and 
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The maximum effective temperature T max = 0.488T, 
occurs at x = 1.36. 

The construction of a disk spectrum proceeds in 
four stages. First, the disk is divided into a set of 
concentric rings or annuli, with each annulus behav- 
ing as an independent plane-parallel radiating slab. 
The vertical structure of each ring is computed, using 
the program TLUSDISK (Hubeny 1990a, b), which 
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is a derivative of the stellar atmosphere program 
TLUSTY (Hubeny 1988). The computation treats 
the entire vertical structure of each ring, from the 
midplane to the surface, without any artificial intro- 
duction of the notion of a distinct "disk atmosphere" . 
The disk stratification is computed in the LTE ap- 
proximation. Hydrostatic equilibrium is assumed, 
with a height-varying (tidal) gravitational accelera- 
tion. Energy balance is enforced between radiative 
losses at the disk surface and heat generation due to 
viscosity, the latter distributed throughout the verti- 
cal extent of the disk. It is assumed that there is no 
flux incident on the disk surface from outside. 

In the specific calculations presented in this pa- 
per, the number of depths between disk mid-plane 
and disk surface is 99. The viscosity prescription is 
based on a Reynolds number approach (Lynden-Bcll 
& Pringle 1974; KHz & Hubeny 1986). A value of 
Re = 5000 is assumed for all radii. A parameter £ in- 
troduced by KHz & Hubeny (1986) controls the verti- 
cal distribution of viscous heat input. In the present 
work, ( — 2/3. The contribution of microturbulence 
to the gas pressure is assumed to be zero in the present 
work. 

Second, the spectrum synthesis program SYN- 
SPEC described by Hubeny, Lanz, & Jeffery (1994) is 
used to solve the radiative transfer equation to com- 
pute the local, rest-frame spectrum for each ring of 
the disk. In addition to detailed profiles of the H 
and He lines, the spectrum synthesis includes "metal" 
lines up to Z = 28 (Ni), dynamically selected from the 
extensive line list of Kurucz & Bell (1995; see also 
Kurucz 1991). Typically, many thousands of lines are 
included in the rest-frame spectrum of each ring. At 
this stage, the wavelength array is irregular and dif- 
fers for each ring. In addition to the angle-averaged 
Eddington flux H(\), specific intensities I(X;fi) at a 
number of emergent angles i are computed and stored 
in this step. 

Third, the rest-frame intensities are combined to 
generate an integrated disk spectrum, using the pro- 
gram DISKSYN5a as described more fully in §3. 
Fourth, the monochromatic fluxes f\ are convolved 
with an adopted Gaussian instrumental broadening 
function and resampled uniformly in wavelength. 

2.2. Choice of Grid Parameters 

In our framework, a steady-state disk model is 
specified by choice of the parameters M wd , R w d, m, 
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Re, and C- These are supplemented by data that de- 
termine the details of how the programs TLUSDISK 
and SYNSPEC compute the atmosphere and spec- 
tra. This section describes the selection of the main 
parameters of the grid. 

Disk models have been computed on a grid of 
M w( j and m with almost-regular spacing but irreg- 
ular grid boundaries. The choice of M w( j is assumed 
to specify i? w d through the mass-radius relation for 
cold electron-degenerate matter. We use an analyti- 
cal mass-radius relation from P. P. Eggleton (private 
communication; cf. Nauenberg 1972). Table [l] lists for 
each choice of M w( j the corresponding white dwarf ra- 
dius i? w( j, surface gravity log<? (in cgs units), Keple- 
rian orbital speed at the surface v SUI { , and the combi- 
nation M w( j/i?^, d which enters into the formula for 
given above. The values of m were chosen to increase 
by factors of VTO 3.16, and the intent was to have 
Af w d / R^d increase by the same ratio from one choice 
of M w d to the next, so that the same T* would apply 
to several different combinations of M w( j and rh in the 
grid; this was achieved for all except the last step in 
Af w d, where the factor is only 3.01. 

The original grid extended only to rh < 1CP 9 Mq yr 
but was later extended. Thus the labeling of models 
by letters, shown in Table [| is slightly irregular. This 
table gives for each disk model: the white dwarf mass 
Af w d, the mass transfer rate m, the number of annuli 
computed, the maximum effective temperature of the 
disk (occurring at x = r/R w d — 1-36), the Rosseland 
optical depth r c to the disk midplane at that disk ra- 
dius, the effective temperature and optical depth for 
the outermost (last) ring, and the radius of the disk. 

For each disk model, stratified atmospheres and 
spectra were computed at a set of radii that increase 
approximately logarithmically from the innermost ra- 
dius, which is always x — R m m/Rwd = 1-05. The 
spacing of rings was chosen to ensure a sufficiently 
dense sampling of T e s(r). The outermost ring was 
chosen to that T c ff(-Rout) would be near 10,000 K, 
since cooler zones at larger radii contribute in only a 
minor way to the far- and mid-UV disk flux. (This 
point is discussed further in §4.4 below.) Table|| gives 
the radii and effective temperature structure for a 
number of representative models. For other models 
the dimensionless radii x are the same. The tem- 
perature structure of other models can be estimated 
closely from the tabulated model that has nearly the 
same T c g at x = 1.36 (see Table |[ or calculated di- 
rectly from the formulas given above. 



The operation of program TLUSDISK is controlled 
by a data file that, among other things, specifies 
the details of opacity and equation of state calcula- 
tions. In our models H and He always contributed 
to the equation of state; in some models cooler than 
12,000 K, the first 30 elements of the periodic ta- 
ble were used. Hydrogen bound-free continua aris- 
ing from the eight lowest energy levels were always 
included in the opacity calculation. Continua from 
fourteen levels of He I were included for models hotter 
than 12,000 K, and eight levels of He II were also in- 
cluded for models hotter than 25,000 K. Other opacity 
sources included electron scattering and free-free con- 
tinua, and the H~ bound-free and free-free continua 
when applicable. For models cooler than 16,000 K, 
the atmosphere computation was carried out with 53 
frequency points defining the continuum; for models 
between 16,000 K and 30,000 K, 59 frequency points 
were used; for models with T c g > 30, 000 K, 65 fre- 
quency points were employed. 

For each ring of each disk model, the atmosphere 
calculation was started with a gray LTE disk model 
(Hubeny 1990a) in hydrostatic equilibrium, then it- 
erated to thermal equilibrium. The convergence cri- 
terion was that the last relative change in any model 
quantity be smaller than 10~ 3 . 

3. SPECTRUM SYNTHESIS FOR RINGS 
AND DISKS 

The radiative transfer solution to compute the de- 
tailed emergent spectrum from each ring was done 
by SYNSPEC at essentially full resolution. That is, 
the maximum spacing between adjacent wavelengths, 
which is 0.01 A for A < 1300 A and 0.02 A for 
A > 1300 A, is nearly small enough to resolve the 
thermal Dopplcr width of metal lines. Microturbulcnt 
broadening of lines was assumed to be zero. The spac- 
ing of computed wavelengths in SYNSPEC is done to 
ensure that there is a wavelength point at the cen- 
ter of each line that contributes to the spectrum, and 
at least one wavelength point between each adjacent 
pair of lines. Thus the wavelength grid is irregular at 
this stage of computation. 

For each ring, both angle-averaged fluxes H v and 
specific intensities l v are computed. The l v compu- 
tation is done for ix = cosz = 1.00,0.75,0.50, and 
0.25. I u is converted to I\ in subsequent computa- 
tional steps. 

To compute the integrated disk spectrum, a sum- 
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mation over rings is carried out. Each ring is divided 
into a large number of azimuthal sectors, and the con- 
tributions are summed with appropriate area weight- 
ing and Doppler shifting. The size of the Doppler 
shift (due to projected orbital motion of the gas) and 
the specific intensity both vary depending on the in- 
clination i of the disk symmetry axis to the line of 
sight. For arbitrary values of fi = cosi, the intensity 
I\ (A; /i) is interpolated linearly in /i between the pre- 
computed I\(X;fj,). Here I\(X;fi = 0) is assumed to 
be zero. 

The resulting disk spectrum is sampled rather 
finely, with a uniform sampling interval of 0.05 A. 
Since the orbital Doppler motions of the gas in the 
disk have smeared each line substantially for all but 
exactly face-on orientations, this high sampling fre- 
quency is redundant. The spectra are convolved with 
an instrumental profile and sampled with four wave- 
length points per FWHM. Except when resolution is 
specifically discussed, all of the spectra described here 
or available as machine-readable files have been con- 
volved with a Gaussian profile whose FWHM is 1.0 

A. 

Disk spectra have been computed and tabulated 
for six inclination angles i. Table ^ provides a handy 
correspondence between /z, i, and sini, which is the 
factor by which orbital motions in the disk are pro- 
jected on the line of sight. The six angles are sufficient 
to show the dependence of the observed spectrum on 
the inclination, but may not always be sufficient to 
allow direct interpolation for the purposes of detailed 
fits of data. Denser sampling of i or computation 
of synthetic spectra for specific "off-grid" angles is 
straightforward; interested readers may contact the 
first author. 

The disk spectra are presented for a "non-projected" 
disk viewed from a distance of d = 100 pc. The cal- 
culation thus amounts to approximating the flux in- 
tegral 

dr 



2k 



/a(A;/x) = 



x / dXh(X,r,iJ,)V(r,<f>,n)A 
Jo 



' 1 f GM wd \ 1/2 

1 H — I sin i sin < 

c V r 



Here 5 is the Dirac delta- function, so that the 
Doppler-shifting function A selects the appropriate 




rest-frame wavelength A' that contributes at the ob- 
server's wavelength A. Therefore the spectra do cor- 
rectly take account of projected velocities. 

The visibility factor V(r, <j>, /i) in the formula makes 
allowance for the presence of an opaque central body 
(the white dwarf) that occults a portion of the inner 
disk. That is, some of the innermost sectors do not 
contribute to the integrated disk spectrum, if they 
would be blocked from view by the white dwarf. This 
visibility factor is not the same as the visibility factor 
for a "flared" disk, as discussed by Robinson et al. 
(1995) or Robinson, Wood, & Wade (1998). In the 
present work flaring effects are neglected (but see the 
discussion in §4.1). 

It is important to note that the results illustrated 
in this paper, as well as the tabulated spectra, do not 
include a factor cos i representing the geometric fore- 
shortening of the (flat) disk, that would normally ap- 
pear as a multiplicative factor in the expression above. 
This factor is omitted, so that the full effect of limb 
darkening alone can be judged from the figures. 

To be compared with actual data, the present 
synthetic disk spectra must be supplemented with a 
model spectrum of the central object, itself with a 
suitably disk-occulted lower hemisphere. This white 
dwarf spectrum is omitted from the results shown 
here. The size of the white dwarf's contribution to the 
total system light is discussed briefly in §4.3. Nguyen 
et al. (1998) have presented some examples of syn- 
thetic solar-abundance white dwarf spectra for the 
mid-UV. A large grid of theoretical solar-composition 
white dwarf model atmospheres has been constructed, 
and will be published in due time (Hubeny & Lanz, 
in prep.). The interested reader may contact I.H. for 
more information and selected models. 

ASCII tables of the synthetic disk spectra are avail- 
able in machine-readable form by anonymous ftp from 
ftp.astro.psu.edu in directory pub/wade/disks. 
There are 26 different files, each corresponding to one 
disk model (as listed in Table ||) . Each file contains a 
single header line, followed by 4601 data lines. Each 
data line contains seven columns: wavelength fol- 
lowed by non-projected flux fx in erg cm~ 2 s^ 1 A -1 
for six values of /i. The initial and final wavelengths 
are 850.00 A and 2000.00 A, and the wavelength spac- 
ing is 0.250 A. Table || shows excerpts from the file 
zz . grid for disk model zz. Publications making use 
of these tables should make bibliographic citation of 
the present paper. 
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4. SOME PROPERTIES OF THE RING AND 
DISK SPECTRA 

4.1. Disk Vertical Structure 

A detailed discussion of the stratification of disk 
models is left to another paper. A couple of issues 
deserve mention, since they relate to the reliability of 
the approximations made in computing the present 
spectra. 

All the rings in all the models have Rosscland 
optical depths r c ;> 1, where r c is measured verti- 
cally, from the outside of the disk to the disk mid- 
plane. (Note that the Rosseland depth referred to is 
computed using only the continuous opacity sources.) 
This high optical depth helps to ensure that: (1) 
the local photosphere, defined as a geometrical region 
where most of the observed spectrum features orig- 
inate, i.e., where the monochromatic optical depths 
measured from the surface are of the order of unity, is 
a relatively thin layer on top of the bulk of the disk, 
consistent with the plane-parallel approximation; (2) 
the local atmosphere can be considered as being lat- 
erally homogeneous, also consistent with the plane- 
parallel assumption; and (3) the formation of (weak) 
absorption lines is not unduly affected by shear due 
to differential rotation of the disk, consistent with 
the hydrostatic approximation. In disk models with 
lower values of m than considered here, some rings 
may become optically thin in the continuum, so that 
the thermal balance of the bulk of the disk is dom- 
inated by line cooling. In that case, the details of 
radiative transfer in the lines, such as turbulent or 
shear broadening, become overwhelmingly important, 
and the strong lines are in emission. In the present 
models, it is expected that such effects arc at least 
not dominant for the computed absorption spectrum. 
The lowest values of the midplane optical depth are 
attained in the inner disk, inside the maximum of 
T e g. They are lowest for models with the lowest m 
at a fixed M wd , or for the highest M wd at a fixed 
to, reaching r c ps 1 for m = 10 -10 ' 5 M yr _1 and 
M wd = 1.210 M© (model x) at x = 1.05. 

The assumption that the photosphere of the disk 
is planar was made in combining the ring spectra to 
make the integrated disk spectra. As will be discussed 
elsewhere, the opening angles of the disk are typically 
9 ~ h/r < 0.1 radian. For the most extreme in- 
clinations discussed here, differential visibility of the 
"front" and "back" halves of the disk, combined with 
a non-negligible difference in limb darkening of the 



front and back halves, should be taken into account 
(see, e.g., Robinson et al. 1995; Wade 1996; Robinson, 
Wood, & Wade 1988). However, this is most impor- 
tant in considering situations in which the disk is ef- 
fectively resolved, such as by eclipse mapping, rather 
than for the integrated light of disks. Another cau- 
tion is that the present models do not include any 
radiation from an outer disk "rim" . 

4.2. Properties of Spectra of Individual Rings 

Because the disk atmospheres considered here are 
optically thick, there is a fairly close resemblance be- 
tween the emitted rest-frame spectrum from some 
disk annulus and the spectrum from a (plane-parallel) 
stellar atmosphere, having the same T e g and with 
log g corresponding to the effective gravity at the pho- 
tosphere of the disk. Small details will be different, 
since hydrostatic equilibrium is solved in one case un- 
der conditions of constant gravitational acceleration, 
and in the other case with a gravitational acceleration 
that increases approximately linearly with distance 
from the midplane. Another potential difference be- 
tween the disk and atmospheric structure arises be- 
cause the total flux is constant with depth in the case 
of stellar atmospheres, but departs from constancy 
in disks. On the other hand, both stellar and disk 
atmosphere models solved in LTE have similar tem- 
perature stratification: in both cases the local tem- 
perature monotonically decreases outward, so major 
features of the spectrum, such as absorption line for- 
mation and limb darkening, are shared in common. 

(To be more general, the local temperature de- 
creases vertically in disks, only when the viscosity 
power-law exponent £ is chosen to be larger than 
about 0.5 for typical cases. For lower values, i.e., 
when the kinematic viscosity does not decrease suf- 
ficiently rapidly towards the disk surface, one faces 
the so-called disk thermal catastrophe [e.g. Shaviv & 
Wehrse 1986]. This is a sharp increase of local tem- 
perature in the superficial layers, where energy is still 
generated by viscous dissipation, while all radiative 
cooling transitions are optically thin and therefore in- 
effective in carrying away the generated energy. With 
our choice of C this catastrophe does not occur.) 

Figure [l] illustrates the computed rest-frame spec- 
tra of three rings (18, 22, 26) from model jj, for the 
region 1500 - 1600 A. The effective temperatures of 
these rings are 29150 K, 19650 K, and 13140 K. Angle- 
averaged Eddington fluxes, \ogH\, are shown. The 
left panel (a) shows the spectra at the full computed 
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resolution, while the right panel (b) shows the same 
spectra convolved with a Gaussian instrumental pro- 
file (FWHM = 0.2 A, corresponding to 39 km s" 1 ) 
and sampled every 0.05 A. Even with this very modest 
instrumental broadening, it is clear that information 
is lost. In a disk spectrum, however, the situation is 
even worse: the Keplerian speed of these rings would 
be ~ 1565, 1180, and 890 km s^ 1 , respectively, so the 
blending of lines due to a mixture of projection fac- 
tors v(r) sin i sin <j> is quite severe, even at low inclina- 
tion. For this reason the disk spectra discussed below 
are convolved with a wider Gaussian (FWHM = 1.0 
A) without significant loss from the point of view of 
comparison with observation. Blending is discussed 
further, below. 

Figure |l| also illustrates the dependence of certain 
spectral features on effective temperature. Notable 
are the Si II lines near 1530 A, the C IV doublet at 
1550 A, and the group of lines near 1560 A (mainly 
C I, Si I, and Fe II), each of which forms at a char- 
acteristic range of temperatures. In a disk spectrum, 
each of these spectroscopic features will be kinemat- 
ically broadened by an amount characteristic of its 
radius in the disk. The integrated disk spectrum thus 
encodes, in a complicated way, the T e g(r) structure of 
the disk. Decoding this scrambled message necessar- 
ily requires spectrum synthesis of the kind presented 
here. 

4.3. Properties of Integrated Disk Spectra 

Figure || shows the far-UV spectrum of a nearly 
face-on disk, model bb, viewed at inclination i = 8?1 
(/i = 0.990), for three different instrumental resolu- 
tions (Gaussian FWHM = 0.1, 1.0, and 3.5 A). The 
smallest FWHM corresponds to a resolution slightly 
worse than that of the ORFEUS Berkeley spectrom- 
eter or the FUSE spectrometers, while the largest 
FWHM corresponds to the HUT spectrometer. For 
this extreme case of minimal Doppler broadening 
and blending of lines, it is clear that improving the 
FWHM from 3.5 A to 1.0 A results m an increase 
of usable information, while further improving the 
FWHM to 0.1 A does not result in significant new 
information about line formation in the disk. Indeed 
the difference between the 0.1 A and the 1.0 A cases 
is barely noticeable at the scale of the figure, in the 
cores of the lines. This shows again that the adopted 
FWHM = 1.0 A is a reasonable choice. 

At a fixed distance from the observer, a disk can 
appear brighter if it is made hotter at constant size 



(higher T max from higher fa) , if it is made larger while 
keeping T max the same (lower M w( j, hence larger i? w( j, 
requiring an increase in m), or if it is tilted closer to 
a face-on orientation (increased projected area, and 
diminished limb darkening) . A disk will generally ap- 
pear "bluer" (relatively more flux at shorter wave- 
lengths) if T max is higher, or if \x is larger (limb dark- 
ening is stronger at shorter wavelengths) . If Tmax and 
/it are kept constant along a disk model sequence of 
changing M w< a (such as the sequence x, t, p, k, aa), 
the effects on the spectrum are expected to be more 
subtle, due to different effective gravities and different 
orbital speeds (hence different line blending) for cor- 
responding rings in different models. Some of these 
effects are illustrated below, while others are better 
studied with the aid of the machine-readable files. 

Figure ^ shows, for a fixed viewing angle /i = 0.75, 
how the flux at 1455 A varies with changes in the 
mass transfer rate, along sequences of constant white 
dwarf mass and radius, and along sequences of con- 
stant T max (recall T max = T e g at x = 1.36). At 
constant M w( j, increasing rh raises the effective tem- 
perature of each disk ring, and also raises the fluxes 
from the disk as a whole. At fixed T max , decreasing 
M w d increases the flux from the disk, since the scal- 
ing of temperature with dimensionless radius x is not 
changed, but the linear scale of the disk increases with 

-Rwd- 

It was noted earlier that the model spectra pre- 
sented here do not include the white dwarf's contribu- 
tion to the light (but do account for occultation of the 
inner disk by the white dwarf). As a rough guide for 
comparing the white dwarf flux with the disk fluxes 
shown in Figure ^|, an (unocculted) white dwarf of ra- 
dius 10 9 cm, viewed from a distance d = 100 pc, pro- 
duces a flux at 1455 A of log f v « 0.5, 0.8, 1.0, 1.1, and 
1.2 (mJy) for effective temperatures of 25000, 30000, 
35000, 40000, and 45000 K, respectively. Fluxes for 
other white dwarf radii can be derived by appropri- 
ate scaling. For an opaque disk with i? m i n = i? w d, the 
projected geometrical area of the white dwarf is re- 
duced by a factor (l + cosz)/2 owing to occultation by 
the disk. The white dwarf may be competitive with 
the disk in the mid-UV, provided it is sufficiently hot, 
and either M w( j or M is sufficiently small, or if the 
disk is viewed at a high inclination angle. 

Figure |J shows, also for fi = 0.750, the behavior 
of two disk flux ratios, /„(1075 A)//„(1455 A) and 
/„(1455 A)// 1/ (1945 A), as m is varied. The former 
ratio relates far-UV and mid-UV flux, while the latter 
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ratio is a measure of the mid-UV slope of the spec- 
trum. As before, solid lines connect models at the 
same M w( j, while dashed lines connect models hav- 
ing the same r max - At T max increases, the flux ra- 
tios (or colors) become "bluer". The mid-UV slope 
is relatively insensitive to disk properties, and above 
Tmax ~ 39, 000 K shows hardly any change. The far- 
UV/mid-UV ratio is more sensitive to T max , as ex- 
pected. Note that the flux ratios are sensitive to disk 
inclination, through the effect of limb darkening, dis- 
cussed below. 

The previous two figures illustrated the behavior 
of disk fluxes for a fixed viewing direction, /i = 0.75. 
Even apart from effects resulting from the changing 
projected area, the flux from a disk will change with fi 
owing to limb darkening. The colors (i.e., flux ratios) 
will change as well, since limb darkening is wavelength 
dependent. These effects are illustrated in Figure [|, 
which shows log /„ vs [i for two disks, models bb and 
z, at three wavelengths f075 A, f455 A, and f945 A. 

Limb darkening occurs for both stars and disks. 
Invoking the classical Eddington-Barbier relation, the 
specific intensity ~ S v {r v — jS), where S v is the 

source function, given in LTE by the Planck function, 
B V (T). Since the temperature decreases (vertically) 
outward, the specific intensities at lower values of fj, 
sample lower-temperature regions, thereby producing 
a spectral energy distribution that is both redder and 
dimmer . In the case of disks there is no averaging 
over a hemisphere as there is in the case of stars, 
because disks lack spherical symmetry. The full effect 
of limb darkening is therefore observable, and disks 
appear very considerably dimmer and reddened when 
viewed at high inclination. Diaz, Wade, & Hubeny 
(1996) have discussed ultraviolet limb darkening for 
disks in detail. 

Attention is now turned from photometric prop- 
erties to the appearance of the actual disk spectra. 
Figure ^ shows a sequence of far-UV fx spectra for 
disks around a M wc i = 0.80 M Q white dwarf, viewed 
with fi — 0.50, normalized near f 330 A. The quantity 
that varies is the mass transfer rate m. The dramatic 
change in overall spectral slope is due to the increas- 
ing temperature of the disk. The hydrogen lines, too, 
show changes in strength and profile. Two effects are 
at work for the lines: (f ) as the disk becomes hotter 
overall, relatively larger areas of the disk are formed 
at high temperatures where the H lines are weak; (2) 
as the disk is made hotter, the location in the disk 
at which the H lines attain maximum strength moves 



outwards, toward lesser orbital speeds and smaller as- 
sociated Doppler broadening. 

The change in line shape due to the change in 
Doppler broadening is also noticeable in the weaker 
lines. For example, near 1135 A note how "absorp- 
tion" metamorphoses into "emission" at the highest 
m. This is actually a case of separate features com- 
bining in the less luminous models to form a single 
absorption feature, placed roughly halfway between 
their rest wavelengths, but only if the disk is cool 
enough overall that these features are formed in the 
inner disk where orbital speeds are high. (See also 
Figure 

Figure |7] shows the mid-UV spectra for the same 
sequence of disk models, again viewed at fi = 0.50. 
The normalization of spectra is now at 1950 A. Hotter 
disks are "bluer" (more flux at short wavelengths), 
and the detailed line spectra again show sensitivity 
to (1) the effective temperatures attained in the disk 
and (2) the velocities characterizing the radii in the 
disk where those temperatures are attained. 

Another view of the far-UV spectrum is provided 
in Figure |[ where a single disk model, bb, is viewed 
from a variety of different lines of sight labeled by /it. 
Recall that the geometric projection factor cos i = fj, 
is not included, so the change in flux levels shown 
in the figure is due to limb darkening alone. Limb 
darkening is very strong in the ultraviolet, and is 
sensitive to both temperature and wavelength. Also 
evident is the strong dependence of the shape and 
strength of features on inclination angle. Since dif- 
ferent fi's correspond to different velocity projection 
factors, the pattern of minima and maxima in the 
disk integrated spectrum changes markedly from one 
example to the next. Individual features appear to 
broaden, split, merge, or shift their locations. The 
overall spectrum becomes smoother when the disk is 
viewed more nearly edge-on. 

Note in Figure || how local maxima in the disk 
flux may appear at the rest wavelengths of strong fea- 
tures, when the feature has been split in two by the 
Doppler effect, such as the 1085 A line of He II, viewed 
at fj, = 0.500. Lines formed at high temperature in 
the inner disk will be broadened more, at a particu- 
lar inclination, than lines formed in the more slowly 
moving, cool outer disk. Thus a single "broadening 
convolution" cannot be applied to a disk spectrum, as 
it can be to the spectrum of a single star: spectrum 
synthesis is mandatory. 
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A final exploration of the dependence of disk spec- 
tra on temperature is provided by Figure ^, which 
shows the region 1300 - 1700 A for a sequence of disks 
with increasing rh onto a 0.80 Mq white dwarf. The 
viewing angle is nearly face-on in this diagram, and 
fluxes are plotted logarithmically to facilitate compar- 
ison of feature strengths in one model with another. 
"Cooler" disks (models m, n, p) show certain absorp- 
tion features more strongly than hot disks, and vice 
versa. Thus, the integrated disk spectra do retain, in 
some measure, information concerning the local, rest- 
frame spectra (cf. Figure [j]) of the rings from which 
they are constituted. This information appears in a 
diluted and convoluted form, which makes interpre- 
tation of the spectra difficult. However, it is the most 
direct connection to the physics of a disk atmosphere. 

An important question is whether the ultravio- 
let spectrum of a steady-state accretion disk is suf- 
ficiently distinctive to enable the white dwarf mass, 
mass transfer rate, and inclination to be recovered 
uniquely. The preceding discussion shows that in 
many cases it will be possible to make useful distinc- 
tions among the models. However, it is possible to 
find pairs of spectra that are very similar. Figure [l^ 
compares two models, cc and v which have the same 
value of T max , achieved through different combina- 
tions of M w( j and rh. Since M w( j differs, there is a 
difference in the Kcplcrian orbital speed of gas in the 
inner disk, but different viewing angles can be chosen 
to match the projected velocities of the disk gas, about 
3400 km s _1 in this case. Likewise, the fluxes from 
the two model disks can be scaled to match at some 
wavelength, as if choosing the distances from which 
the disks are viewed. The result of this matching, as 
shown by the Figure, is two spectra that are nearly 
indistinguishable. The remaining differences arc due 
to different amounts of "reddening" from limb dark- 
ening and a small difference in the effective photo- 
spheric gravity of the two disks. The differences are 
subtle and might pass unnoticed or be attributed to 
differences in interstellar reddening, for example. Yet 
the difference in rh is more than a factor of three! 

Additional, rather precise, knowledge about an ob- 
served system may be required to break such a de- 
generacy of model spectra that match the spectro- 
scopic observations about equally well. This informa- 
tion might concern the orbital inclination, the white 
dwarf mass, or the distance. The specific match cho- 
sen for illustration does not involve matching extreme 
values of mass, inclination, or distance, however, and 



if encountered in practice would require precision of 
measurement better than is usually attained, in order 
to distinguish the models. 

4.4. Missing Flux Due to Truncation at T(r) » 

10,000 K 

An important consideration in using these models 
is to what extent flux is missing, due to the fact that 
the spectrum synthesis was stopped when T e g fell be- 
low roughly 10,000 K. This missing flux may on the 
one hand affect the overall "color" of the integrated 
spectra, and may on the other hand affect the details 
of line blending and feature strength. The expecta- 
tion is that, although the outer disk (i.e., r > -Rout of 
the present models) is large in surface area, it is cool 
enough that the contribution in the far- and mid-UV 
spectral region can be neglected. 

The missing flux question was investigated quan- 
titatively for seven models (m, s, t, u, v, dd ; jj) 
representing the whole range of T max of the disk 
model grid. From visual inspection of all the indi- 
vidual ring spectra of model jj, three wavelength in- 
tervals were selected, each 10 A wide and centered 
at 1075 A, 1455 A, and 1945 A. In these intervals, 
the line blocking is relatively small and either consis- 
tent or smoothly varying over the range of temper- 
atures that was modeled. For each spectral interval 
and each of the seven disk models, mean fluxes H v 
and intensities l v were computed, ring by ring. Plots 
of x 2 H u {x) and x 2 I v {x) against In a; (x = r/i? w( j) 
were constructed, and these were used to extrapolate 
x 2 H u {x) or x 2 I v (x) smoothly to zero in the outer disk. 
Numerical integrations of iftotai = 271"/ xH v {x)dx = 
2tt J x 2 H l/ (x)d(lnx) over the extrapolated disk and 
over the actually computed disk limits were then com- 
pared. A similar procedure was used with I v . 

In this way it was estimated that at 1075 A, an 
error is made that is never larger than 0.2 per cent 
for any disk model, either for angle-averaged flux or 
for specific intensity. For 1455 A, the maximum error 
grows to 3.5 per cent (for the coolest overall model, 
model m). At 1945 A, the maximum estimated miss- 
ing flux is largest, as would be expected, amounting 
to 9.3 per cent for model m, but no larger than 3.1 per 
cent for models as hot or hotter than model u. The 
error in intensity is largest when the disk is viewed 
nearly face-on (jj, = 1). At large inclination angles 
limb darkening, which is stronger at a fixed wave- 
length for cooler atmospheres, diminishes the relative 
importance of the missing outer disk. 
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As a further check, angle-averaged fluxes H\ for 
solar-abundance stellar atmospheres (Kurucz 1994) 
were used to synthesize approximate disk spectra cor- 
responding to models m, s, t, u, v, dd and jj. These 
were truncated at T out ~ 10000 K (as with the self- 
consistent disk models described here), and alterna- 
tively at T out « 4000 K. Comparison of the disk 
fluxes gives truncation errors directly for these stel- 
lar atmosphere-based model disk spectra. These er- 
rors are very similar to the truncation errors for the 
self-consistently computed disk spectra of this paper, 
as estimated by the extrapolation method described 
above. 

5. DISCUSSION 

Modern observations of CVs in the ultraviolet are 
frequently of good enough quality that they deserve 
to be compared with model spectra that are more 
detailed than models used in the past. The grid of 
model spectra presented here is intended to provide a 
basis for such comparisons. The construction of the 
spectra involves many decisions, such as what level of 
approximation is to be used, or what values or ranges 
of parameters are to be considered, so these topics 
have been addressed at some length. Also, the illus- 
trative comparison of models has been emphasized, 
to show the dependence of the spectrum on m, /i, etc. 

Detailed analysis of observed CV spectra is beyond 
the scope of this presentation, as it involves a host of 
additional considerations pertaining to individual ob- 
jects. These may include such things as interstellar 
reddening and the extent to which the white dwarf 
mass or orbital inclination may be constrained from 
orbital studies. It is appropriate, however, to review 
briefly some issues that may arise during the com- 
parison of model spectra and observed spectra. We 
emphasize that the models presented here, although 
advanced by comparison with past efforts, and cer- 
tainly suitable for use, also represent a starting point 
for improved future models. 

The present disk spectra do not include the contri- 
bution from the central white dwarf or its boundary 
layer. The boundary layer is expected to contribute 
mainly in the extreme ultraviolet, A < 911 A (Wade 
1991; Polidan, Mauche, & Wade 1990). The white 
dwarf's contribution to the far- and mid-UV spec- 
trum can be estimated (see §4.3), and can be readily 
computed and included. 

The present disk models and spectra are based on 



vertical hydrostatic equilibrium. Observations show 
that among luminous CVs, many of the resonance 
lines and stronger subordinate lines in the UV exhibit 
P-Cygni or blue-shifted absorption profiles, indicating 
that a wind exists. Likely, this wind arises from the 
disk, rather than from the central star (Pereyra, Kali- 
man, & Blondin 1997; Proga, Stone, & Drew 1998). 
The direct observational consequence of the wind is 
that several of the more obvious spectral diagnostics 
cannot be used, that would otherwise provide infor- 
mation about the properties of the disk. More sub- 
tle problems may also arise, if the weak UV lines are 
formed in gas that already participates in the outflow, 
or if the spectrum is formed in an extended zone that 
undermines the plane-parallel (i.e., laterally homoge- 
neous) approximation that was used in computing the 
models of disk annuli. Further modeling will be re- 
quired to address these questions. 

The disk stratifications computed here, and the 
spectra derived from them, do not include the effect of 
irradiation by an external source, whether it is direct 
illumination of the disk by the central star, illumina- 
tion of the outer disk by the hotter inner disk, or emis- 
sion or scattering from a wind. Such irradiation will 
be most important for relatively cool zones of the disk 
(Hubeny 1990b), for which the contrast between the 
color temperature of incident radiation and the pho- 
tospheric gas temperature is largest. Since the far- 
and mid-UV spectra of disks are dominated by the 
inner disk, it is reasonable as a first step to neglect 
this complication, provided that the disk is viewed 
as a whole. If spectra from isolated portions of the 
disk are considered, irradiation becomes a larger is- 
sue, even in the UV. 

Finally, within the class of hydrostatic, non-irradiated 
models there are some important questions to be 
further explored, concerning, e.g., line blanketing, 
or non-LTE treatment of stratification and radiative 
transfer. These issues are common to all computa- 
tions of astrophysical atmospheres, whether stars or 
disks. The model spectra presented here serve as 
a reference, against which to compare spectra com- 
puted under different assumptions or approximations. 

6. SUMMARY 

Based on the need for detailed models to com- 
pare with modern observations of disks in CVs, we 
have presented a grid of far- and mid-UV spectra 
of steady-state, hydrostatic, non-irradiated accretion 
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disks. The grid covers a significant range of M w( j 
and to, suitable for comparison with novalike vari- 
ables or dwarf novae in outburst. The grid consists 
of 26 model disks, with spectra computed for each of 
them at six inclination angles. The models were com- 
puted using computer programs TLUSDISK, SYN- 
SPEC, and DISKSYN. A basic set of continuum opac- 
ities was used in the calculation of vertical stratifica- 
tion, but a large line list was used in the computation 
of the emergent spectra. The emergent spectra take 
account of limb darkening, and also of Dopplcr shifts 
due to orbital motion of the gas in the disk around the 
central white dwarf. Occultation of part of the inner 
disk by the central star is also taken into account. 

The spectra extend from 850 A to 2000 A, suitable 
for comparison with IUE, HST, HUT, ORFEUS, and 
FUSE spectra. The models extend to cool enough 
temperatures that essentially all of the light produced 
by the disk at the wavelengths of interest is accounted 
for. Machine-readable files accessible via anonymous 
ftp provide tabulations of the spectra at 0.25 A steps 
with a FWHM resolution of 1.0 A. The spectra are 
presented in the form of non-projected fluxes for a 
disk viewed from a distance of 100 pc. For compar- 
ison with observations, the disk fluxes may need to 
be supplemented with fluxes from the appropriately 
occulted central star. 

The radiative properties of both individual annuli 
and full disks arc discussed with illustrative examples, 
to show how the disk spectra incorporate information 
about projected gas velocities and the variation of 
T e g with radial distance from the white dwarf. In 
the disks considered, which are optically thick at all 
radii, the local UV spectra resemble spectra of stars 
that share the same T e fj and log g. Fluxes and col- 
ors (flux ratios) of the integrated disk spectra depend 
on M w( j and to. Limb darkening is shown to have 
a strong effect in the ultraviolet, affecting both the 
overall flux level and colors of the disk. Line broad- 
ening and blending is complex and non-intuitive in 
these disks. Synthetic spectra such as those presented 
here may be needed to understand the details of ob- 
served spectra, or even in some cases to identify the 
lines responsible for a given spectral feature. Despite 
trends in the disk spectra that are evident with vary- 
ing accretion rate or inclination, spectra from rather 
different models may nevertheless appear quite simi- 
lar, and additional information about the CV may be 
required to decide on a unique, best model. 

Several possible improvements or additions to the 



present set of models are discussed. As such mod- 
els become available it will be of interest to compare 
the resulting spectra with those presented here. In 
the meantime, the present grid of spectra represents 
a self-consistent treatment of both the vertical strat- 
ification and emergent spectral properties of steady- 
state disks, with sufficient detail in the computation 
of the line spectrum to be of use in analyzing observed 
spectra of CV disks. 

Support from NASA grants NAG5-1698, NAG5- 
2125, NAGW-3171, and NAG5-3459 is gratefully ac- 
knowledged. 
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Fig. 1. — Rest-frame spectra (surface Eddington 
flux H\) for three rings of disk model jj. These 
rings, Numbers 18, 22, and 26, have T off = 29150 K, 
19650 K, and 13140 K respectively. Left panel: spec- 
tra are shown at full computed resolution, maximum 
stepsize = 0.02 A. Right panel: spectra are shown af- 
ter convolution with a Gaussian instrumental profile, 
FWHM = 0.2 A. Different spectral features form at 
different characteristic temperatures. 

Fig. 2. — Non-projected far-UV fluxes from disk 
model bb, viewed at 100 pc distance. The disk is 
viewed nearly face-on. Three instrumental resolutions 
are shown, specified by Gaussian FWHM = 0.1, 1.0, 
and 3.5 A. FWHM = 3.5 A roughly corresponds to 
the resolution of the HUT, while FUSE is expected 
to resolve better than 0.1 A. The disk models do not 
include wind lines. 

Fig. 3. — Non-projected flux f v , averaged over the 
band 1450 — 1460 A, in mJy (One milliJansky = 
10~ 26 erg cm -2 s _1 Hz -1 ). The fluxes refer to disks 
viewed at a distance of 100 pc. Limb darkening and 
line blending appropriate to fx — cosi = 0.750 have 
been taken into account. Points corresponding to disk 
models d, m, x, ff, and jj are labeled. Solid lines 
connect models with the same M w( j (rightmost curve 
0.35 Mq, leftmost curve 1.21 Mq). Dashed lines con- 
nect models with roughly the same T max (lowest curve 
- 16, 500 K, uppermost curve - 92, 000 K). 

Fig. 4. — Ratios of flux /„, for disks viewed from 
an inclination i = 41?4 (fi = 0.750). Points cor- 
responding to disk models d, m, x, and ff are la- 
beled. Solid lines connect models with the same M w <j 
(rightmost curve 0.35 Mq, leftmost curve 1.21 Mq). 
Dashed lines connect models with roughly the same 
T max (lowest curve ~ 16, 500 K, uppermost curve 
~ 92,000 K. The color range (at a fixed value of //) 
is narrower for the flux ratio involving longer wave- 
lengths, as expected. 

Fig. 5. — Solid circles: the non-projected fluxes 
/„ at 1075, 1455, and 1945 A (top to bottom, at 
H = 1) for disk model bb (M wd = 0.550 Mq, rh = 
10~ 8 ' 5 Mq yr _1 ). Open circles: the same for disk 
model z (M wd = 1.210 M ,m = lO" 90 M Q yr" 1 ). 
Model bb has T max = 39110 K, while model z has 
Tmax = 51520 K. Colors (flux ratios) and fluxes 
clearly depend on \i, in addition to M wc i and m. 
Non-projected flux shows the effect of limb darken- 
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ing only; multiply by a factor = cosi to account for 
geometric foreshortening. One milliJansky (mJy) = 
10~ 26 erg cm -2 s" 1 Hz -1 . 

Fig. 6. — Far-UV fluxes fx, normalized at 1330 A, for 
a sequence of disk models [m, n, p, q, cc, hh) with 
constant white dwarf mass but increasing m (bottom 
to top) . All spectra are shown as seen from a viewing 
direction i = 60° (/j, — cosi = 0.50). Spectra are con- 
volved with a Gaussian instrumental profile, FWHM 
= 1.0 A. 

Fig. 7.— Mid-UV fluxes fx, normalized at 1950 A, 
for the same sequence of disk models shown in Fig. 6. 
Details are the same as for Figure 6. 

Fig. 8. — Far-UV, non-projected fluxes fx from disk 
model bb, (M wd = 0.550 M©, m = 10~ 8 - 5 M© yr^ 1 ), 
for six inclinations, \l — 0.99, 0.95, 0.75, 0.50, 0.25, and 
0.15 (top to bottom). The /i = 0.99 curve is offset ver- 
tically by 0.5 x 10~ n erg cm -2 s _1 A^ 1 and shown 
with a heavier line for clarity. Effects of limb dark- 
ening and line blending are apparent. Fluxes should 
be multiplied by a factor // = cos i to account for 
geometrical foreshortening. 

Fig. 9. — Details of the mid-UV spectra of a sequence 
of disk models, showing the characteristic features of 
disks as a function of disk temperature T max . The 
spectra are shown as seen from a viewing angle of 
i = 8?1 so that line blending is minimized. Non- 
projected fluxes are for a viewing distance of 100 pc. 
As with all figures in this report, the flux from the 
white dwarf is not included. 

Fig. 10. — A comparison of two closely matched disk 
spectra. Each model disk (v, cc) has the same run 
of T G ff with dimensionless radius cc, and the viewing 
angles fi = cosi have been chosen to match approx- 
imately the run of projected velocities v(x)sini. Fi- 
nally, the fluxes have been normalized at 1950 A (up- 
per panel) or 1330 A (lower panel). The spectrum 
of model v is plotted with a heavier line for clarity. 
Residual differences between the spectra are due to 
limb darkening and a systematic difference in effec- 
tive gravities amounting to about 0.17 dex (see Table 
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Table 1 

Adopted White Dwarf Parameters 



M wd 


Rwd 


log g 


Vsurf 


M wd /B? wd a 


(M e ) 


(10 9 cm) 




(km s _1 ) 




0.35 


1.142 


7.55 


2020 


0.2374 


0.55 


0.905 


7.95 


2840 


0.7423 


0.80 


0.699 


8.34 


3890 


2.347 


1.03 


0.518 


8.71 


5140 


7.423 


1.21 


0.378 


9.05 


6520 


22.35 



a Units are M /(1O 9 cm) 3 . 
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Table 2 
Parameters of Disk Models 



(EXTERNAL PLANOTABLE, PORTRAIT MODE) 



Table 3 

Temperature Structure of Selected Disk Models 



(EXTERNAL PLANOTABLE, LANDSCAPE MODE) 
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Table 4 

Inclination Angles for Computed Disk Spectra 



smt 



/i = COS I 



0.990 
0.950 
0.750 
0.500 
0.250 
0.150 



8?1 
18?2 
41?4 
60? 
75? 5 
81?4 



0.141 
0.312 
0.661 
0.866 
0.968 
0.989 



Table 5 
Sample Table for Model zz 



(EXTERNAL PLANOTABLE, LANDSCAPE MODE) 
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1000 1050 1100 1150 

Wavelength (A) 



TABLE 3 

Temperature Structure of Selected Disk Models 



Ring 


x = r/R w d 


Model m 


Model n 


Model p 


Model q 


Model cc 


Model hh 


Model jj 






Teff 


T e ff 


T e ff 


T e ff 


T e ff 


T e ff 


T e ff 


01 


1.05 


12840 


17130 


22840 


30450 


40610 


54160 


72220 


no 


i on 


IDUzU 


zloOU 


oc/1 nn 


qcnnn 
ooUUU 


OUDDU 


O / OOU 


nnnnn 

yuuyu 


03 


1.36 


16490 


22000 


29330 


39110 


52160 


69560 


92760 


04 


1.56 


16180 


21580 


28780 


38380 


51180 


68250 


91010 


05 


1.80 


15450 


20610 


27480 


36650 


48870 


65170 


86910 


06 


2.05 


14620 


19500 


26010 


34680 


46240 


61670 


82240 


07 


2.35 


13680 


18240 


24320 


32440 


43250 


57680 


76920 


08 


2.70 


12700 


16930 


22580 


30110 


40150 


53450 


71400 


09 


3.10 


11730 


15650 


20860 


27820 


37100 


49480 


65980 


10 


3.55 


10819 


14430 


19240 


25660 


34210 


45630 


60840 


11 


4.10 


9900 


13200 


17600 


23470 


31300 


41740 


55660 


1 


A Id 
4. (\j 




ioi nn 


i £1 /in 


01 f.on 


osfinn 




t^i nQn 


13 


5.40 




11060 


14750 


19660 


26220 


34970 


46630 


14 


6.20 




10090 


13460 


17950 


23930 


31910 


42560 


15 


7.10 






12290 


16390 


21850 


29140 


38860 


16 


8.20 






11140 


14860 


19820 


26430 


35240 


17 


9.40 






10140 


13530 


18040 


24060 


32090 


18 


10.80 








12290 


16390 


21860 


29150 


19 


12.40 








11160 


14880 


19850 


26470 


20 


14.30 








10100 


13460 


17960 


23940 


21 


16.40 










12220 


16300 


21730 


22 


18.90 










11050 


14730 


19650 


23 


21.70 










10010 


13350 


17800 


24 


25.00 












12060 


16080 


25 


28.80 












10890 


14530 


26 


33.10 














13140 


27 


38.00 














11880 


28 


43.70 














10740 



1 



3 




1.5 

1 .8 .6 .4 .2 

/u ( = cos i) 



TABLE 5 

Excerpt from ASCII Table zz. table for Model zz 



wavelength 




mu=0.150 




mu=0.250 


mu=0.500 




mu=0.750 




mu=0.950 


mu=0.990 


850.000 


2 


.21100E-12 


3 


.78900E-12 


4.68300E-12 


5 


.39800E-12 


5 


.84900E-12 


5.91300E-12 


850.250 


2 


.21100E-12 


3. 


.78700E-12 


4.68100E-12 


5 


.39600E-12 


5 


.84600E-12 


5.90100E-12 


850.500 


2 


21000E-12 


3. 


.78500E-12 


4.67800E-12 


5 


.39300E-12 


5 


.84200E-12 


5.88700E-12 


850.750 


2 


.20900E-12 


3. 


.78300E-12 


4.67500E-12 


5 


.38900E-12 


5 


.83600E-12 


5.87500E-12 


851.000 


2 


.20800E-12 


3 


.78100E-12 


4.67200E-12 


5 


.38500E-12 


5 


.82900E-12 


5.86800E-12 


851.250 


2 


.20700E-12 


3 


.78000E-12 


4.67000E-12 


5 


.38200E-12 


5 


.82100E-12 


5.86500E-12 


1998.750 


1 


.06600E-12 


1 


.81100E-12 


2.33700E-12 


2 


.74100E-12 


3 


.00400E-12 


3.05800E-12 


1999.000 


1 


.06600E-12 


1 


.81000E-12 


2.33700E-12 


2 


.74100E-12 


3 


.00800E-12 


3.06400E-12 


1999.250 


1 


.06600E-12 


1 


81000E-12 


2.33600E-12 


2 


.74000E-12 


3 


01200E-12 


3.06700E-12 


1999.500 


1 


.06500E-12 


1 


.80900E-12 


2.33500E-12 


2 


.74000E-12 


3 


.01600E-12 


3.06900E-12 


1999.750 


1 


.06500E-12 


1 


80900E-12 


2.33400E-12 


2 


.73900E-12 


3 


.01900E-12 


3.06900E-12 


2000.000 


1 


.06400E-12 


1 


.80900E-12 


2.33300E-12 


2 


.73900E-12 


3 


.02100E-12 


3.06800E-12 




1400 1600 1800 2000 

Wavelength (A) 
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0.80 M©, /u = 0.990 (i = 8.1°) 




log dM/dt = -10.5 -> -8.5 



14 i i i i i i i i i i i i i i i i i i i i i 

1300 1400 1500 1600 1700 

Wavelength (A) 







_ 1 1 1 1 1 


1 1 


I I I I _ 




Model 


v, /u = 0.75 — 


— 


Model 


cc, ju = 0.50 

— 


1 1 1 1 1 1 


i i 


i 1 i i i 



1.5 — 



.5 




900 



1400 1600 1800 

Wavelength (A) 



2000 



Model v, ii = 0.75 

Model cc, yu = 0.50 




1000 1100 1200 

Wavelength (A) 



1300 



TABLE 2 

Parameters of Disk Models 





M. .j 

lv± wa 


log 771 


Rings 


rp a 

± max 


T b 

' C 


iu & y 


rp d 

± out 


iu & y 


loff T\ .. . + 

±U B ±L out 




\ 1V1 Q) 


\ ivi y j 














(pin 1 


d 


0.350 


-9.5 


i-ii 


16490 


93 


5.72 


9900 


4.86 


9.67 


e 


0.350 


-9.0 


1-14 


22000 


200 


5.85 


10090 


4.70 


9.85 


aa 


0.350 


-8.5 


1-17 


29330 


216 


5.98 


10150 


4.52 


10.03 


n 


U.ooU 


— o.U 


i on 




con 


a no 


i ni on 
1U1UU 




lU.zl 


h 


0.550 


-10.0 


i-ii 


16490 


46 


5.85 


9900 


4.95 


9.57 


j 


0.550 


-9.5 


1-14 


22000 


76 


6.00 


10090 


4.84 


9.75 


k 


0.550 


-9.0 


1-17 


29330 


90 


6.14 


10150 


4.67 


9.93 


bb 


0.550 


-8.5 


1-20 


39110 


178 


6.26 


10100 


4.47 


10.11 


gg 


0.550 


-8.0 


1-23 


52160 


350 


6.38 


10100 


4.27 


10.29 


m 


0.800 


-10.5 


1-11 


16490 


8 


5.95 


9900 


5.02 


9.46 


n 


0.800 


-10.0 


1-14 


22000 


22 


6.12 


10090 


4.95 


9.64 


P 


0.800 


-9.5 


1-17 


29330 


41 


6.29 


10150 


4.80 


9.82 


q 


0.800 


-9.0 


1-20 


39110 


65 


6.43 


10100 


4.62 


10.000 


cc 


0.800 


-8.5 


1-23 


52160 


138 


6.55 


10100 


4.43 


10.18 


hh 


0.800 


-8.0 


1-25 


69560 


230 


6.67 


10890 


4.36 


10.30 


s 


1.030 


-10.5 


1-14 


22000 


7 


6.21 


10090 


5.03 


9.51 


t 


1.030 


-10.0 


1-17 


29330 


20 


6.42 


10150 


4.92 


9.69 


u 


1.030 


-9.5 


1-20 


39110 


28 


6.59 


10100 


4.77 


9.87 


V 


1.030 


-9.0 


1-23 


52160 


54 


6.72 


10100 


4.58 


10.05 


dd 


1.030 


-8.5 


1-25 


69560 


105 


6.84 


10890 


4.53 


10.17 


jj 


1.030 


-8.0 


1-28 


92760 


160 


6.97 


10740 


4.31 


10.36 


X 


1.210 


-10.5 


1-17 


28970 


7 


6.48 


10020 


4.99 


9.55 


y 


1.210 


-10.0 


1-20 


38640 


14 


6.72 


9970 


4.86 


9.73 


z 


1.210 


-9.5 


1-23 


51520 


22 


6.88 


9890 


4.70 


9.91 


zz 


1.210 


-9.0 


1-25 


68710 


48 


7.00 


10760 


4.68 


10.04 


ee 


1.210 


-8.5 


1-28 


91620 


79 


7.13 


10610 


4.47 


10.22 



a Effective temperature at x = 1.36. 

b Midplane Rosseland mean optical depth r, at x = 1.36. 
c logg measured at r = 1, x = 1.36. 
d Effective temperature at the outermost ring. 
e logg measured at r = 1, outermost ring. 
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